EXTENSION THEORY FOR 
ELLIPTIC PARTIAL DIFFERENTIAL OPERATORS 
WITH PSEUDODIFFERENTIAL METHODS 



GERD GRUBB 

Abstract. This is a short survey on the connection between general extension 
theories and the study of realizations of elliptic operators A on smooth domains 
in K", n > 2. The theory of pseudodifferential boundary problems has turned 
out to be very useful here, not only as a formulational framework, but also 
for the solution of specific questions. We recall some elements of that theory, 
and show its application in several cases (including new results), namely to 
the lower boundedness question, and the question of spectral asymptotics for 
differences between resolvents. 



1. Introduction 

The general theory of extensions of a symmetric operator (or a dual pair of 
operators) in a Hilbert space, originating in the mid-1900's, has been applied in 
numerous works to ordinary differential equations (ODE), and also in a (smaller) 
number of works to partial differential equations (PDE). 

There is a marked difference between the two cases: In ODE, the playground 
for boundary conditions is usually finite-dimensional vector spaces, where linear 
conditions can be expressed by the help of matrices. Moreover, the domains of 
differential operators defined by closure in i2-based Hilbert spaces can usually all 
be expressed in terms of functions with the relevant number of absolutely continuous 
derivatives. 

In contrast, boundary conditions for PDE (in space dimensions n > 2) are pre- 
scribed on infinite-dimensional vector spaces. Moreover, the domains of differential 
operators in L2-based spaces will contain functions with distribution derivatives, 
not continuous and possibly highly irregular. 

Whereas extensions of ODEs can usually be described in terms of matrices, the 
tools to interpret extensions in PDE cases are therefore much more complicated. We 
shall give a survey of some tools developed through the years, and their applications, 
emphasizing the use of pseudodifferential operators. 

Outline. In Section 2 we recall the basic issues of elliptic boundary value problems. 
Pseudodifferential operators (?/>do's) are introduced in Section 3, and Section 4 in- 
troduces pseudodifferential boundary operators (i/>dbo's). In Section 5 we recall the 
elements of a general abstract extension theory, and in Section 6 we show how this 
is implemented for realizations A of an elliptic operator A on a domain SI C K™. 
Section 7 focuses on the resolvent formulas that can be obtained via the general 
theory. In the last sections we go through several cases where pseudodifferential 
techniques have proved extremely useful (some of the results here are quite recent) : 
In Section 8 it is the question of whether lower boundedness holds simultaneously 
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for a realization A and the operator L over the boundary that enters in the cor- 
responding boundary condition; the new results deal with unbounded domains. In 
Section 9 it is the question of showing Weyl-type spectral asymptotics formulas for 
differences between resolvents. The results there go back to the early theory, and 
Section 10 presents some additional new results. 

2. Elliptic boundary value problems 

In the following we use the customary multi-index notation for differential op- 
erators: d = d x = (di, . . . , d n ), dj = d Xj = d/dxj, and D = D x = (Di, D n ), 
Dj = D Xj = -id/dxf, then d a = d" 1 ■ ■ ■ d%", D a = I?" 1 • • • D%», for a G Nft; here 
|a| = ai H h a n . 

A differential operator of order m > 0, 

A= a a (x)D a 

\a\<m 

is said to be elliptic, resp. strongly elliptic, on a set U C W l , when the principal 
symbol 

a m (x,0 = Y a <*( x )C 

\a\—m 

satisfies 

a m {x,0 ^ 0, resp. Rea m (x,^) > 0, 

for x G U , £ G W 1 \ {0}. As a basic example, the Laplacian A = d\ + ■ ■ ■ + d% has 
principal symbol (and symbol) equal to — 1£| 2 , so it is elliptic, and —A is strongly 
elliptic. The Laplacian has been studied for several hundred years, and the problems 
around it solved by explicit solution formulas. It is the cases with variable (in- 
dependent) coefficients, and domains more general than simple geometric figures, 
that have been a challenge in more modern times. 
The problem 

(1) An = f 

for a given function / on a subset of R™ usually has infinitely many solutions. 
To get a problem with unique solvability, we must adjoin extra conditions such 
as suitable boundary conditions. We can consider A with a domain consisting of 
the functions satisfying the boundary condition, as an operator A acting between 
suitable spaces. 

Then the question of existence of a solution corresponds to the question of 
whether A is surjective, and the question of uniqueness of a solution corresponds 
to the question of whether A is infective. In this way, the question of solvability of 
differential equations is turned into a question of properties of specific operators. 
The operator point of view became particularly fruitful when it was combined with 
appropriate scales of function spaces, such as the Sobolev spaces, Sobolev [So50], 
and with Distribution theory, Schwartz [Sch50]. 

When 17 is a smooth open subset of K™ with boundary dfl — E, we refer to the 
standard Z/2-Sobolev spaces, with the following notation: H s (M. n ) (s G R) has the 
norm \\v\\ 8 = WJ 7 ^ 1 ((^} s J r v)\\ L2 ^y, here T is the Fourier transform 

T: u(x) ^ (J"m)(0 = = / e- ix <u{x) dx, 
and (£) = (1 + |£| 2 )^- Next, with m denoting restriction to O, 



H s (n) =r Q H s (R n ), 
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provided with the norm ||u|| s = inf{||v|| s | v G H s (W l ), u = rnv}. Moreover, 

iJ s (H) = {ue H s (R n ) | suppu C O}, 

closed subspace of H s (R n ); it identifies with the antidual space of H~ s (Cl) (the 
space of antilinear, i.e., conjugate linear, functionals) , with a duality consistent 
with the L 2 scalar product. For s equal to a nonnegative integer k, Hq(Q) is 
usually written Hq(Q). Spaces over the boundary, i/ s (E), are defined by local 
coordinates from H s (M."^ 1 ), s G M. (There are many equally justified equivalent 
choices of norms there; one can choose a particular norm when convenient.) When 
s > 0, there are dense continuous embeddings 

ff'(£) C L 2 (E) C H- S (Z), 

and there is an identification of H~ S (I]) with the antidual space of H S (T,), such that 
the duality (ip, i(j)- StS coincides with the L2(E)-scalar product when the elements 
lie there. Detailed explanations are found in many books, e.g. Lions and Magcnes 
[LM68], Hormander [H63], [G09]. (There is a difference of notation: For k G N, 

the spaces denoted H 2 (f2) in [LM68] are not the same as our H 2 (SI) that are 
consistent with [H63]; they have the best duality and interpolation properties.) 

Consider the case where A is defined on a smooth open subset SI of R™ and 
has coefficients in C°°(Q), and assume that A is elliptic on CI. The results in this 
case are a model for results under weaker smoothness hypotheses. One defines the 
maximal realization A max as the operator acting like A in the distribution sense 
with domain 

D(A max ) = {ue L 2 (Q) I Au G L 2 (fl)}- 
it is a closed, unbounded operator in L 2 (S1). The minimal realization A m j n is defined 
as the closure of A acting on C£°(Cl) (the compactly supported C°°-functions on 
Cl). When Vl is bounded, or is unbounded and there are suitable bounds on the 
coefficients of A, 

D(A mhl ) = H™(fl). 
The formal adjoint A' of A is the differential operator acting as follows: 

A'u= D a (a a (x)u). 

\a\<m 

By definition, A' min and A max are adjoints of one another (as unbounded operators 
in L 2 (fl)). 

The linear operators A satisfying 

G ^4 G ^4max 

are called realizations of A. 

Generally, A max is far from being injective, whereas A m - m is far from being sur- 
jective, but it may be possible to find realizations A that are bijective from D(A) 
to L 2 {0). 

We see that the theory of distributions (which in this context was preceded by 
the definition of differential operators acting in the weak sense) allows defining 
operators representing the action of A in a generalized sense. Here invertibility can 
sometimes be achieved by methods of functional analysis. A fundamental example 
is the Dirichlet problem (where jqu = u|s) 

(2) Au — f in Vl, 70M = <p on E, 

for a strongly elliptic second-order operator having Ke(Av, v) > c||v||^ 2 ^ with 
c > for v G Cq°(CI). By use of the so-called Lax-Milgram lemma one could 
define a realization A 1 of A with D{A 1 ) C Hq(Q), such that A 1 : D(A 1 ) — > L 2 (fl) 
bijectively. (Details are found in many books, e.g. [G09], Ch. 12.) 
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But then the question was: How close is A^ 1 to solving the problem in a more 
classical sense? Second-order derivatives have a meaning on H 2 (Q), by closure of 
the definition on C 2 (il), so one can ask: 

• if / e L 2 (n), is u e h 2 (Q)? 

• More generally, if / e H k (Ct) for some k e N , is u e H k+2 (9)1 

The answer was first found for the behavior of u in the interior of O: Indeed, 
when / G H k (fl), u is in i/ fe+2 over subsets of O with positive distance from the 
boundary. This is the so-called interior regularity. 

There remained the question of regularity at the boundary. It was answered 
positively in papers by Nirenberg [N55] and Ladyzhenskaya (see the account in 
[L85]). This was followed up by research on higher-order operators A and more 
general boundary conditions Tu = ip (possibly vector valued), where results on 
interior regularity and regularity at the boundary were established under suitable 
conditions. Besides ellipticity of the operator A one needs a condition on how the 
boundary condition fits together with A. Some authors called it the "covering con- 
dition" or the "complementing condition" , but the name "the Shapiro-Lopatinskii 
condition" (after [Sh53], [L53]) has been more generally used. It is also customary 
to call the system {^4, T} elliptic when it holds (this was suggested in Hormander 
[H63]; we return to a motivation in Section 4). A fundamental reference in this con- 
nection is the paper of Agmon, Doughs and Nirenberg [ADN59] that collects and 
expands the knowledge on elliptic boundary value problems. An important point 
of view was to obtain so-called "a priori estimates" (estimates of a Sobolev norm 
on u by norms on Au and Tu plus a lower order norm on it), shown for smooth 
functions at first, and extended to the considered solution. 

Important monographs exposing the theories and the various authors' own con- 
tributions were written by Agmon [A65], Lions and Magenes [LM68]; the latter 
moreover contains valuable information on the surrounding literature. The early 
theory is exposed in Courant and Hilbcrt I— II [CH53], [CH62]. 

In the works at that time, although the striving to show existence of a solution 
operator was always in the picture, the emphasis was more on showing qualitative 
properties of the unknown function u in terms of properties of the given data / and 
ip, regardless of whether u could be described by an operator acting on {/, ip} or 
not. 

Direct machinery to construct approximate solution operators in general came 
into the picture with the advent of pseudodifferential methods. 

3. Pseudodifferential operators 

One of the few cases where an elliptic differential operator has an explicit solution 
operator is the case of I — A on R", whose action can be described by use of the 
Fourier transform T as (1 — A)u = J 7-1 ((1 + £ | 2 ) Tu) , and whose solution operator 
is 

°p(nV)«=^'(TTWH 

A variable-coefficient elliptic differential operator on R" can also be described by 
the help of the Fourier transform, 

Au = ^2 a a (x)D a u= ^2 ftaW^^r^") = ^^(s^)^ = Op(a(i,?))?j, 

|a|<m |a|<m 

where a(x, £) = X^| a |<m a a( x )£" i s the symbol. But even when the symbol satisfies 
a ( x >0 f° r au the operator 

Op^er 1 ) = F^aix^)- 1 ? 
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is not an exact inverse. Nevertheless, it is useful in the discussion of solutions, since 
one can show that 



where the remainder 1Z is of order —1 (lifts the exponent of a Sobolev space by 1). 

A thorough treatment of Op(a(s,() _1 ) and suitable generalizations that come 
closer to being an inverse of A (such "almost-inverses" are called parametrices) 
was obtained with the invention of pseudodifferential operators, -0do's. Some of the 
initiators were Mihlin [M48], Calderon and Zygmund [CZ57], Seclcy [Sc65], Kohn 
and Nirenberg [KN65], Hormandcr [H65], [H71]. 

General f/'do's are defined from general symbols p(x, £) as 



where we use the notation d~£ — (2ir)~ n dt;; here p(x,£) is required to belong to a 
suitable class of functions. 

Not only do the operators make sense on K" where the Fourier transform acts, 
they are also given a meaning on manifolds, by use of coordinate change formulas 
and cutoff functions. The theory is not altogether easy; it uses concepts from 
distribution theory in a refined way. Moreover, it is not exact but qualitative in 
many statements, so it can be something of a challenge to derive good results from 
its use. A fine achievement is that it leads to Fredholm operators, when applied to 
elliptic operators on compact manifolds. Here one is just a small step away from 
having invertible operators; this can sometimes be achieved by relying on additional 
knowledge of the situation. 

A so-called "classical" ipdo is an operator defined from a symbol that has an 
asymptotic series expansion in homogeneous terms (a polyhomogeneous symbol): 



P(X,0 ~ Pm-j(x,0^ Pm-j(x,t£) = t m 3 p m -j(x,0 for |£| >1,*>1. 



It is said to be of order m, and Op(p) maps H s to H s ~ m for all s£l. p m is called 
the principal symbol, and p is said to be elliptic when p m {x,£) f° r l£l > 1- 
Here one has that 



where Hi and IZ2 are of order m + m' — 1. In this way, the principal part dominates 
the behavior. When p is elliptic, the principal part of a parametrix is found as p" 1 
(for |£| > 1, extended smoothly to £ e R n ). Also the notation p° is used for the 
principal symbol. 

There exist many different symbol classes with generalizations of the above prop- 
erties, designed for particular purposes. 

The very attractive feature of classical -0do's is that they form a scale of operators 
of all integer orders, including differential operators among those of positive order, 
and including parametrices and inverses of elliptic differential operators among 
those of negative order. Moreover, it is an "algebra" , in the sense that the elements 
by composition (and by addition) lead to other classical ^do's. 

The calculus is explained in the original papers and in several subsequent books, 
such as Treves [T80], Hormander [H85]; a detailed introduction can be also found 
in Chapters 7-8 of [G09]. 

The tpdo theory gives (after one has done the work to set it up) an easy proof 
of interior regularity of solutions to elliptic problems. 



Op(a( a ;,£))Op(a(x,e)- 1 ) =1 + 11, 




o P (p) op(p') = o P (pp') + Ki = Op( PmP ' m ,) + n 2 
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4. PSEUDODIFFERENTIAL BOUNDARY OPERATORS 

When an elliptic differential operator A is considered on a subset of R n or on 
a manifold with boundary — let us here for simplicity just consider the case of 
a smooth bounded open subset of W 1 — we must impose boundary conditions 
to get uniquely solvable problems. Let assume that we are in a case where the 
boundary condition Tu = cp together with (1) gives a uniquely solvable problem; 
here T is a trace operator mapping functions on il into M-tuples of functions on 
E = d£l. We can formulate this in terms of matrices: 

/a\ c°°(n) c°°(n) 

(3) :C°°(fi)^ x has an inverse (R K) : x -> C*°°(ft). 

\t J c°°(e) m c°°(e) m 

Here K is called a Poisson operator; it solves the semi-homogeneous problem 

Av = in O, Tu = 95 on E. 
The operator R solves the other semi-homogeneous problem 

Aw = f in SI, Tw = on E. 
In a closer analysis of i?, we can write it as a sum of two terms: 

(4) R = Q+ + G, 

where Q is the ipdo A -1 on R™, Q + is its truncation r + Qe + to 0, and G is a sup- 
plementing operator adapted to the specific boundary condition, called a singular 
Green operator (s.g.o.). The operator e + stands for extension by (to functions on 
R"), and the operator r+ stands for restriction to O. 

The calculus of pseudodifferential boundary operators (V'dbo's) was initated by 
Boutet de Monvel [B71], who introduced operator systems encompassing both the 
systems (^) and their solution operators (i? K). The original presentation is 
somewhat brief, and was followed up by extended expositions, in the detailed book 
of Rempel and Schulze [RS82] , which elaborated the index theory, and in the paper 
[G84] which completed some proofs of composition rules (with new points of view), 
and showed spectral asymptotic estimates for singuar Green operators. The book 
[G96], whose first edition was issued in 1986, developed a calculus of parameter- 
dependent V'dbo's, leading to resolvent and heat operator constructions. The recent 
book [G09] gives a full introduction to the theory. 

In the systematic calculus of Boutet de Monvel one consider systems (called 
Green operators): 

/P+ + G K\ C°°(f7) w C 00 ^)^' 
A = : x -> x , where 

\ T S) C°°(E) M C°°(E) M ' 

• P is a ps.d.o. on W 1 , satisfying the so-called transmission condition at E 
(always true for operators stemming from elliptic differential operators); 

• P + = r + Pe + is the truncation to f2 (the transmission condition assures 
that P+ maps C°°(Ti) into C°°(fi)); 

• T is a trace operator from f2 to E, K is a Poisson operator from E to O, S 
is a V'do on E; 

• G is a singular Green operator, e.g. of the type KT. 

The composition of two such systems is again a system belonging to the calculus. 

The operators extend to act on Sobolev spaces. For T and G there is a condition 
expressing which differential trace operators -fjU = (9/9n) J u|s that enter: T or G 
is said to be of class r when only 7/s with j < r enter; and then they act on H S (Q) 
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for s > r — \ . The class case is the case where they are purely integral operators, 
well-defined on L<2,(Si). 

All entries can be matrix-formed. They are defined in local coordinates by 
formulas involving Fourier transformation and polyhomogeneous symbols. The idea 
is as follows: In local coordinates at the boundary, where fi and £ are replaced by 
K+ = {x G R" | x n > 0} and M™" 1 (with points x' = (x 1: . . . ,x n _i)), the system 
has for each (x',£') a boundary symbol operator acting in the a;„-variable: 

a(x',Z',D n ) = 

(5) /p(x',0,?,D n )+g(x',li',D n ) k(x',?,D n )\ H m (R+) N L 2 (R+) N ' 

V t(x',?,D n ) s(x',?) ) ' C M C M ' 

Here m is the order of the operator. Each entry in a acts in a specific way. E.g., 
when the matrix is (^) in (3), the boundary symbol operator is the model operator 
obtained by freezing the coefficients at x' and replacing derivatives in D", by their 
Fourier transforms (with respect to x 1 G R" -1 ). The principal boundary 

symbol operator a°(x', D n ) is formed of the top order terms. The principal 
boundary symbol operator for (i? is the inverse of the principal boundary 

symbol operator for (^). (For (5), g and t must be of class < m.) 

From the boundary symbol operator one defines a full operator by applying the 
ipdo definition in the x'-variable, 

Op'(a(x',£',D n ))u = J e lx ' e a(x' , , D^iFy'^uiy' , x n )) tig. 

The symbols have asymptotic series of terms that are homogeneous in (£',£ ra ) 
(different rules apply to the different ingredients, and we must refer to the men- 
tioned references for further details). One then defines A to be elliptic, when 

(1) P is elliptic, i.e. its principal symbol p°(x, £) is invertible at each (x, £) with 

Kl > i, 

(2) the principal boundary symbol operator a°(x', D n ) is invertible at each 
{x',£) with \C\ > 1. 

For a system (^) formed of an elliptic differential operator A and a differential trace 
operator T, 2) is precisely the old covering/complementing/Shapiro-Lopatinskii 
condition for {A,T}. 

In the elliptic case, one can construct a parametrix B° from the inverses of the 
symbols in l)-2), such that AB° — I and B°A — I have order < —1, and the con- 
struction can be refined to give errors of arbitrarily low order. With supplementing 
information it can be possible to obtain an inverse. 

For example, there holds a a solvability theorem for an elliptic differential oper- 
ator problem as in (3), formulated in this framework as follows: 

Theorem 4.1. Let il C K" be a smooth, bounded open set, denote d£l — S, and let 

A = J2\ a \<2m a a( x ) Da with a a G C°° (tl) be elliptic on il, i.e., J2\ a \=2 m a a( x )£> a ^ 
for x G O, £ G M. n \ {0}. Let T — (Tj) , JL 1 be a column vector of trace operators 
Tj = foBj, where the Bj are differential operators of order mj with C°° -coefficients, 
< mi < ■ ■ ■ < m m < 2m — 1. (Then T is of class r = m m + 1 < 2m.) Assume 
that {A, T} is elliptic. 

The operator A = (^) defines a continuous mapping 

/A\ H°(Cl) 

(6) A=\ \:H 2m+s (n)^ x fors>r-2m-\, 
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and there is a system B = (R G) (a parametrix) belonging to the calculus and 
continuous in the opposite direction, such that 

H s (n) H s '(n) 

Tlx : x -> x 

TZ 2 : H 2m+S {tt) H 2m+S '(n) 

for all s as in (6), aZZ s' > s. _ffere if is a row vector of Poisson operators (Kj)jL 1 
of orders —mj, and R = Q + + G, where Q is a parametrix of A on a neighborhood 
of fl, and G is a singular Green operator. 

The operator A in (6) is Fredholm for each s, with the same finite dimensional 
kernel and cokernel in C°° for all s. 

If A is bijective, the inverse belongs to the calculus ( and is of the same form as 
B). 

When r = 2m, the lower limit for s is — |; cases where it is < — | occur for 
example for the Dirichlet problem, where r = m, and s can go down to — m — ^. It 
is useful to know that the Poisson operator K in fact has the mapping property 

m 
3 = 1 

for all s G K. The trace operator T is called normal, when joBj = bjj mj + 
J2k< mj Bjk"fk with an invertible coefficient bj for each j. (More general normal 
boundary value problems are described below in Section 9.) 

For example, for a second-order strongly elliptic operator with a Dirichlet con- 
dition, the operator in the theorem maps as follows: 

A=\ : H 2+S (Q) -> x for s > -f, 

with parametrices and solution operators continuous in the opposite direction. 

Elliptic operators A of odd order occur mainly as square matrix-formed oper- 
ators, and there is a similar theorem for such cases, where also the Bj can be 
matrix-formed. Operators of Dirac-type are a prominent first-order example. The 
matrix case is also interesting for even-order operators. The results can moreover 
be worked out for operators defined on manifolds, acting in vector bundles. (See 
e.g. [G74], on the even-order case, for notation and the appropriate definition of 
normal boundary conditions.) 

5. Extension theories 

We shall now recall some elements of the functional analysis theory of extensions 
of given operators. This has a long history, with prominent contributions from J. 
von Neumann in 1929 [N29], K. Friedrichs in 1934 [F34], M. G. Krein in 1947 [K47], 
M. I. Vishik in 1952 [V52], M. S. Birman in 1956 [B56], and others. The present 
author made a number of contributions in 1968-74 [G68]-[G74], completing the pre- 
ceding theories and working out applications to elliptic boundary value problems; 
further developments are found e.g. in [G83], [G84], and in recent works. 

At the same time there was another, separate development of abstract exten- 
sion theories, where the operator concept gradually began to be replaced by the 
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concept of relations. This development has been aimed primarily towards applica- 
tions to ODE, however including operator-valued such equations and Schrodingcr 
operators on R"; keywords in this connection are: boundary triples theory, Weyl- 
Titchmarsh m- functions and Krein resolvent formulas. Cf. e.g. Kocubei [K75], 
Vainerman [V80], Lyantze and Storozh [LS83], Gorbachuk and Gorbachuk [GG91], 
Derkach and Malamud [DM91], Arlinskii [A99], Malamud and Mogilevskii [MM02], 
Briining, Geyler and Pankrashkin [BGP06], and their references. In later years 
there have also been applications to elliptic boundary value problems, cf. e.g. Am- 
rein and Pearson [AP04], N.D. Kopachevskh and S.G. Krcm [KK04], Bchrndt and 
Langer [BL07], Ryzhov [R07], Brown, Marietta, Naboko and Wood [BMNW08], 
Gesztesy and Mitrea [GM08], and their references. 

The connection between the two lines of extension theories has been clarified in 
a recent work of Brown, Grubb and Wood [BGW09]. 

At this point we should also mention the recent efforts for problems on non- 
smooth domains: Posilicano and Raimondi [PR09], Grubb [G08], Gesztesy and 
Mitrea [GM08, GM11], Abels, Grubb and Wood [AGW11]; here [G08, GM11, 
AGW11] use [G68]. 

In the following, we shall use the notation from [G68]-[G74] and [BGW09]. 

Let there be given a pair A min , A' min of closed, densely defined operators in a 
Hilbcrt space H , such that the following holds: 

C (j4 min ) — . A max , ^4 m j n C (j4 m in) • ^max- 

Let M = {A A min cic A max }. Write Au as Au, when Ae M. 

We assume that there is given an operator A^ s M, the reference operator, with 
G g(A 7 ) (the resolvent set); then 

A ■ r A r A A' r A* r A' 

^min -^7 ^- n maxj ^-min ^- 7 ^- max' 

The case where A m i D — A' min and A 7 is selfadjoint, is called the symmetric case. 
Let Z = ker A max , Z' = ker A' max , and define the basic non-orthogonal decom- 
positions 

D(A max ) = D(A 7 )+Z, D{A' m ^) = D(A;)+Z' } 

denoted u = u^ + u^ = pr 7 u + pr^ u, where pr 7 = J 4 7 1 J 4 max , with a similar notation 
with primes. 

By Tpr x u = ux we denote the orthogonal projection from H to X. The injection 
X H is denoted ix (it is the adjoint of pr x : H — > X). 

There holds an "abstract Green's formula" for u € D (A max ), v e D(A' max ): 

(7) (Au,v) - {u,A'v) - ((Au)z>,v c ) - {u c ,{A'v) z ). 
It can be used to show that when A e M, and we define 

V = pi c D(A), W = pr c , D(A*), 

then 

(Au) w } I u € D(A)} is a graph, 
defining an operator T from D{T) C V to W. 

Theorem 5.1. (^G68]J There is a 1-1 correspondence between the closed operators 
A e M and the closed densely defined operators T: V — > W , where V C Z , W C Z' 
( arbitrary closed subspaces ), such that A corresponds to T: V — > W if and only if 

(8) D{A) = {u e D(A max ) | pr c u G D(T), (Au) w = Tpr c «}. 

In this correspondence, V = pr^ D(A), W = pr^, D(A*), and 
• A* corresponds analogously to T* : W — >• V. 
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• keryl = kerT; ranyl = ranT + (H G W). 

• A is bijective if and only if T is so, and then 

A' 1 = A- 1 + iyT- 1 pr w . 

One also has 

D(A) = {u = v + A-\Tz + f) + z\v£ D(A min ), z £ D(T), feZe W}, 
where v, z and f are uniquely determined from u. 

The result builds on the works of Krem [K47] and Birman [B56] (for selfadjoint 
operators) and Vishik [V52], and completes the latter: In Vishik's paper, the A 
were set in relation to operators over the nullspaces going in the opposite direction 
of our T's, and the results were focused on those A's that have closed range (the 
so-called normally solvable realizations). Our analysis covered all closed A. 

The condition in (8) 

(9) (Au) w =TpT ( u 

can be viewed as an "abstract boundary condition" . 

When A £ £?(Ay), one can do the same construction for the operators shifted by 
subtraction of A. We denote 

Z x = kcr(A max - A), Z{ = ker(A^ ax - A), 

and have the decompositions (where pr x = (Ay — \)~ 1 (A max — A)) 

D{A max ) = D{A 1 )+Z X , u = u x + u x =pr x u + pr A u, 

with a similar notation with primes. 

Corollary 5.2. Let A £ g(A^). For the closed A £ M, there is a 1-1 correspon- 
dence 

~ I T x : V\ — > W\ , closed, densely defined 

1 with V\ C Z\, W\ C Z't, closed subspaces. 

Here D(T X ) = pr x D(A), V\ = D{T X ), W x = pr A D(A*), and D(A) consists of 
the functions u £ D(A max ) such that u x £ D(T X ) and 

T x u x = ((A - \)u) W - x . 

Moreover, 

• ker(l - A) = kerT A ; ran(i - A) - ranT A + (H Q W~ x ). _ 

• A — A is bijective if and only if T x is so, and when A £ g(A) n £>(Ay), 

(A A)- 1 = (A, - A)- 1 + iv,^)- 1 pr Wx . 

This gives a Krem- type resolvent formula for any closed A £ M. with q(A) n 
q(A 7 ) ^ 0. _ 

The relation between T and T x was determined in [G74] in the symmetric case, 
for real A, and the proof given there extends immediately to the general situation 
(as shown in [BGW09]): 

For A £ g(A-y), define 

E x =1 + A(Ay - A)" 1 , it has the inverse F x = I - \A~ X , 

and similarly E' x = I + X(A* - A)" 1 has the inverse F' x = I — ^(A*)^ 1 on H. 
Then E X F X = F X E X = I, E' X F' X = F' X E' X = I on H . 
Moreover, E x and E' x restrict to homeomorphisms 

E X :V^V X , E' X :W^W- X , 
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with inverses F v resp. F$. In particular, D(T X ) = E V D(T). 

The operator families derived from E x are related to what was called a gamma- 
field in other works from the 1970's and onwards, as a simple special case. 

Theorem 5.3. Let G vw — — pr w \E x iv ; then 

(E'w)*T x E v = T + G vw . 

In other words, T and T x are related by the commutative diagram 

^A— 3 V 



t+g^ w D(T X ) = E X D(T). 



W x = W 

( E w)* 

In the joint work with Brown and Wood [BGW09] we moreover showed how the 
study relates to studies of boundary triples and M-functions by other researchers 
(as referred to in the start af this section; more references are given in [BGW09]): 

Let V = Z, W = Z', then, with H = Z',IC = Z, the mappings 



Tqu J V U C 



r> 

T' v 



: D(A' max ) -> K x H, 



form a boundary triple: Both mappings ( r *) and ( r , x ) are surjective, their kernels 
are D(A m i n ) resp. D(A' min ), and they satisfy the Green's formula 

(A max u,v) - (u^a^v) = (r lU ,r' v) n - (r u,i»K;, 

which is a rewriting of (7). 

Here one can consider a boundary condition 

(10) rut = TT u, 

where we allow T to be unbounded (closed densely defined) from /C to H; it defines 
a restriction A of A max by D(A) = {u G D(A max ) \ T u G D(T),Tiu = TT u}. 
Then it is customary to define an M-function as follows: 

Definition 5.4. For A G g(A), M(A) : ran(ri— TT ) -> JC is the operator satisfying 

M(A)(ria; - Tr a;) = r a:, for all x G kcr(A max — A) = Z\. 

The analysis in [BGW09] showed that M(X) is a holomorphic family of operators 
in £(H,IC). On the other hand, when A and its boundary condition (10) are 
considered from the point of view of extensions [G68]-[G74] recalled further above, 
A is the operator corresponding to T: Z — > Z' by Theorem 5.1. Then we find 
moreover, in terms of the A-dcpendent families introduced in that context: 

M(A) = -(T + G^z,)- 1 = -F X {T X )-\F' X )*, when A G g(A)n Q (A J ). 

This gives the Krein resolvent formula in the form 

(A- A)- 1 = (A, A)- 1 - i Zx E x M(X)(E' x r pr z , . 

For the case of general V and W, we could likewise construct an M-function 
from W to V for A G g(A), and establish a Krein resolvent formula around it. The 
following result is shown in [BGW09]: 
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Theorem 5.5. Let A correspond to T: V —¥ W. For A £ g(A) there is a well- 
defined holomorphic family M(A) £ £(W, V); 

M(A) = P r c (/ - (I - A) _1 (A max - A)) A^W- 

W^en A G £>(A) n g(A-y), then 

M(A) = -(T + G^)- 1 - —Fy(T x )~ 1 (Fw)*, 

and 

(A-X)- 1 = (A 1 - A)- 1 - i Vx E*M(\)(E*r Ww- X ■ 

To have both T x (for A e g(A~ ( )) and M(A) (for A £ £>(A)) available is an 
advantage, since ker(^4 — A) = kcrT A and ran(A — A) = ranT A + (H Q W^) give 
straightforward eigenvalue information at the poles of M(A) in £>(A 7 ). 

Remark 5.6. The name M-function is consistent with the notation in some papers 
that [BGW09] refers to, but possibly diverges from others (one could also use the 
longer name Weyl-Titchmarsh function). There is a recent publication of Malamud 
[M10] that exposes related resolvent formulas on the basis of [MM02]. (Let us 
remark that [M10] seemingly ascribes a hypothesis of separate surjectiveness of r 
and Ti to the presentation in [BGW09]; this is not so.) 

6. Implementation of the abstract set-up for elliptic operators 

We shall now recall how the abstract theory is applied to a concrete choice of 
elliptic operator A. Here A max and A min arc denned as in Section 2; they are 
closed operators in H = L 2 (£l). In [G68] general even-order operators were con- 
sidered, and the reference operator (called A 1 in Section 6) was taken to repre- 
sent a general normal boundary condition. To simplify our explanation, we shall 
here just consider a second-order strongly elliptic operator A and let A-y stand for 
the Dirichlet realization, mentioned after (2). We have by elliptic regularity that 
D(A 7 ) = H 2 (Q) n Hq(Q), and we can assume that (a constant has been added 
to A such that) the lower bound m{A 1 ) is positive. The lower bound m(P) of a 
operator P is defined by 

(11) m(P) = inf{Rc(Pu, u) \ u £ D{P), \\u\\ = 1} > -co; 

when it is finite, P is said to be lower bounded. 

The trace operator 70 defines a continuous mapping H S (Q) — > i/ s ~2(£) for 
s > — i. We shall also need a more advanced fact, namely that, as shown by Lions 
and Magenes (see e.g. [LM68]), 70 extends to a mapping D(A max ) -> J ff-i(S), and 
defines homeomorphisms 

1Z :Z^H-HZ), 7z ,:Z / 4JT-*(£) ) 

where Z and Z' are the nullspaces of A max and A' max (not contained in H s (il) for 
s > 0). The inverse of jz is consistent with the Poisson operator K 1 solving the 
semi-homogeneous Dirichlet problems (2) with / = 0, in the sense that 

Similarly, the inverse 7^, is consistent with the Poisson solution operator K' 7 solv- 
ing the Dirichlet problem for A' with / = 0, K' = iz'lz' ■ Moreover, with A- 
dependence, 

K x = \ Zx lzl K i = izpz>> 
solve the semi-homogeneous Dirichlet problems for A — A, A' — A, when A £ g(A 7 ). 
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For a closed subspace V C Z, let X = -f V C H 2 (S). Here we denote the 
restriction of 70 by -yv', 

(12) lv :V^X, 

with a similar notation for Y = ^qW and A-dependent cases. The map 7y : V — ^ X 
has the adjoint 7^: X* — > V. Here X* denotes the antidual space of X, with a 
duality coinciding with the scalar product in L 2 (£) when applied to elements that 
come from L 2 (E). The duality is written (tp, (p)x*,x- 
We denote 

(13) K ltX = \ vl y l : X -> V C H; 

it is a Poisson operator when X = H~^{Ti). 

Now a given T: V — > is carried over to a closed, densely defined operator 
L : X -> F* by the definition 

i = (7^)*2V, D{L)= lv D{T); 
it is expressed in the diagram 

V X 



W - — * Y* 

There is a similar definition in the A-dependent case. 

Before formulating the results in a theorem, we shall interpret the abstract 
boundary condition (9) defining the realization A, as a concrete condition in terms 
of L. 

A has a Green's formula (for sufficiently smooth u, v) 

(14) (Au,v)n - (u,A'v) n = (^iu,7 w)e - (7ouXu)e, 
where 

vi = S71, = S71 + .4' 70, 
with a nonvanishing smooth function s and a suitable first-order differential oper- 
ator A' on S. 

Let A £ £>(Ay). In addition to the Poisson operators if^ resp. .fT' solving the 
Dirichlet problems for A — A resp. A' — A, we shall need the Dirichlet-to- Neumann 
operators 

P^ Vl =v,K^ P'l v ,=v[K'~\ 

that map the Dirichlet boundary value into the Neumann boundary value for null- 
solutions. By the composition rules for V'dbo's, they are pseudodifferential opera- 
tors of order 1; moreover, it is known that P~ o Vl is elliptic. 

Theorem 6.1. Define the reduced Neumann trace operator T by 

r = 1/1 - P° 0tVl lo = viA- l A m ^: D(A max ) -> iT*(£). 

It is continuous and surjective, and vanishes on Z. With the analogous definition 
for A' one has the reduced Green's formula. 

(Au,v) - (u,A'v) = (Tu,j v)i _i - (7o«,r'w)_i i, 
valid for all u E D(A nlllx ), v £ D(A' max ). In particular, 

(15) (Au, w) = (ru, 7oio)i._i, when u e D(A max ),w 6 Z' . 
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For A G q(A^) we similarly define 

r A = vi - P^lo = vi(Ay - A)" 1 ^^ - A), 

r' A = *i- p^Uio = v'M; - A)- 1 ^ - a), 

continuous and surjective from D(A max ) resp. D(A' max ) to H^(T,); then there holds 
a reduced Green's formula 

(Au,v) - (u,A'v) = (T A u,7o«)i,_i - (lou,T' x v)_^i , 

for u G D(A max ), v G £»(^ ax ). 

Here -D(j4 max ) is provided with the graph-norm. 

Now let A be a closed operator lying between A min and A max , so A £ *M. The 
abstract boundary condition (9) for A may be written: 

(16) (Au,w) = (Tu c ,w), all w £ W. 

The left-hand side equals (Tu,jow)i _i by (15). The right-hand side equals 
(Tu c ,w) = (T-fy 1 ^, 7jyScH = (L7 u,7 w)y.,y, 

by definition of L (it is used that when £ V, = 7y 1 7yu^ = 7y 1 7ou). 
Hence (16) may be written 

(17) (Tu,j w)i_i = (L70M, iow)y,Y, ah w £ W. 

The injection iy : Y — >• if~5(£) has as adjoint the mapping i y : iJ3(S) — »• y* 
that sends a functional ^ on H~i{Yi) over into a functional iyV on ^ by : 

(iyV; l v 9 )y*,r = f° r a U <P £Y. 

Then (17) may be rewritten as 

iyTu — L-f u, 
or, when we use that r = v\ — P° .^70, 

(18) i^ lU =(L + i^P 7 o ^) 7 ow. 

This is the boundary condition derived from (9), when A corresponds to T: V ^ W 
by Theorem 5.1, carried over to L: X — > Y* by (12). 
Then Theorem 5.1 implies: 

Theorem 6.2. There is a 1-1 correspondence between the closed operators A £ M 
and the closed densely defined operators L: X — > Y* , where X and Y are closed 
subspaces of (£), such that A corresponds to L: X — >■ Y* if and only if D(A) 
consists of the functions in D(A max ) for which 

(19) 7o« G D(L), i* Y v lU = (L + iyP^how. 

In this correspondence, X — joD(A), Y = -foD(A*), and 

• A* corresponds analogously to L* : Y — > X* . 

• kcrl = i V 7^ 1 kcri; rani = 7^ rani + (H Q W), cf. (12), (13). 

• A is bijective if and only if L is so, and then 

(20) A- 1 = A- 1 + iy^L- 1 ^ Ww = A- 1 + K^ X L-\K^ Y )*. 
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Theorems 6.1 and 6.2 are from [G68], except that we have modified the notation 
a little. 

In [BGW09], the subspace cases are treated with insertion of an isometry 
A i : i 2 (S) i?~2(E), that allows replacing X and Y by closed subspaces Xi 
and Yi of £2(2), identified with their dual spaces; then iy is replaced by an or- 
thogonal projection pr Yl ■ 

In the case where Y = H~^{E), i.e., W = Z', the map \y is superfluous, and 
the second condition in (19) takes the form 

When also X = we say that A represents a Neumann-type condition 

(21) v x u = Cj a u; here C = L + P° Q Vi on D(L). 

In this case, L can act like a pseudodifferential operator, namely when C (in the 
condition v\u = Cjou) is a differential or pseudodifferential operator. 

Let us consider a slightly different set-up where C is a given first-order differential 
or pseudodifferential operator on E, and we define A by 

(22) D{A) = {ue D(A max ) I v lU = C lo u], 

where 70 and v\ are considered as mappings from D(A max ) to fl-i(E) resp. 
iJ~5(E). We shall discuss the corresponding operator L: X — > Y* . Since 
{ lo u,viu}: H 2 (n) -> ijl(E) x Hi(S), £>(L) = J D(A) D fli(E). Then since 
if 2(E) is dense in iJ _ 5(E), X = i? _ 5(E). By use of Green's formula (14) it is 
checked that the adjoint A* extends the realization of A' with domain consisting 
of the functions v € H 2 (£l) satisfying 

v x v = C 70V, 

so also Y = i/~5(E). Thus we are in the case of Neumann- type boundary condi- 
tions, so by comparison with (6), it is seen that L acts like C — P° Vl . 

The domain D(L) equals {tp G H~i(Y,) \ (C - P% , Vl )<p G H*(E)}; it may not 
be easy to determine more exactly. Note that L is used as a map from its domain 
in (E) to _ff + 2(E), although it acts like a ipdo of order 1. 

One case is clear, though: If C — P® Q Vl is elliptic of order 1, then Lip e H 2 (E) 

implies ip € i?5(E); in this case D(L) = ff5(E). Moreover, D(A) C H 2 {Q). A 
check of the boundary symbol rules shows that this is precisely the case where the 
system {A,ui — C70} is elliptic. Here we have: 

Theorem 6.3. Let C be a first-order differential or pseudodifferential operator on 
E and define the realization A of A by (22). Then if C — P® v is elliptic, the 
operator L: X — > Y* corresponding to A by Theorem 6.2 acts like C — P® aMl and 
has 

X = Y = ff~5(E), £>(£,) = .fff(E). 
Moreover, D(A) C H 2 (il). Related statements hold for the adjoint A*. 

A Robin condition j^u = 670U, with a smooth function b(x) on E, is elliptic, 
since L acts like b — f-? 0l ^ l5 where -P° 0il/1 is elliptic of order 1 and 6 is of order 0. 
In the case of Theorem 6.3, when L is bijective, the formula (20) has the form 

(23) I" 1 =A- 1 +K J L- 1 (K' 7 )*, 

where all ingredients belong to the ^dbo calculus: K 1 is a Poisson operator, L^ 1 is 
a i/'do on E, (K' 7 )* is a trace operator of class 0, and the composition K 1 L~ 1 {K' 1 )* 
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is a singular Green operator (of class 0). A 1 1 is the sum Q + 
ipdo Q on R™ and a singular Green operator, as in (4). 



G of a truncated 



Let us include a few words on higher-order cases: When A is of order 2m, 
there is a Green's formula generalizing (14), where 70 is replaced by an m- vector 
7 = {70i ■ • • , 7m-i}, and v\ and v[ are replaced by m-vectors of trace operators of 
orders m, . . . , 2m — 1, mapping into products of Sobolev spaces of different orders 
over S. One then gets vector versions of the reduced Neumann trace operators F 
and r', with matrix- formed versions of the Dirichlct-to-Neumann pseudodiffercntial 
operators, but the basic mechanisms in the interpretation are the same. There 
are interesting cases of subspaces X, Y of the products of Sobolev spaces over S, 
where cllipticity considerations are relevant. Details are given in [G68]-[G74] and 
[BGW09] . 

7. Resolvent formulas 

When A e g(Aj), there is a similar representation of A — A in terms of a boundary 
condition defined from an operator L x acting over the boundary. Here it is of 
particular interest to find the connection between L and L x , just as we found the 
connection between T and T x . It turns out that the relation between L and L x is 
simpler: They both go from X to Y* , whereas T resp. T x map between different 
spaces due to the shift from Z to Z\. This holds, since 

D(L) = ■y D(A) = l0 D{A - A) = D(L X ), X = TML) = D(L X ), 

with similar statements for D(L*), D((L X )*) and Y. Then we have: 




T+G? 



The horizontal maps compose as 7y A -Ey = 7Vj ilw-J^ 1 i^w)* = (lw)~ 

L x = 1 y 1 {T + G x VM )( 1 * w y 1 . 
In terms of L x , the boundary condition reads (analogously to (18)): 



so 



(24) 



i* Y v lU = (L x + i^P x ) l0 u, l0 u G D(L). 



Since D(A — A) = D{A) is at the same time defined by the boundary condition 
i^is lU = {L+i^P? (o l/i ) 10 uioT 10 u G D(L), we have that L x +i^P x a ^ = L+i^P* a Pi 
on D(L), so 



L x = L + i^{p 



i0 

70,^1 



P x ) on D(L). 



The last formula is convenient, since P® Vl 



P x Vl can be shown to be bounded 



from (E) to (£); hence L x is a perturbation of L by a bounded operator. 

Also the general M-function defined in Section 5 carries over to an M-function 
on the boundary, a holomorphic family of operators Ml (A) E C(Y*,X) defined for 
A G q{A). 

The results are collected in the following theorem: 

Theorem 7.1. Let A correspond to T : V — >• W as in Theorem 5.1, carried over to 
L: X — > Y* as in Theorem 6.2. For A G q{A 1 ) it is also described by the boundary 
condition (24), and there holds: 
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(i) For A G g(A^), £° ^ - P^ Vi e C(H~$ (£), (£)) and 

i A = L + - P^J on D(L). 

(ii) For A e £>(A), £/iere is a related M -function e C(Y*,X), 

M L (X) = 70 (J - (A - A) _1 (^ max - A))A; 1 i w7 ^. 

(iii) For A G £>(A) n g(A 7 ), 

M L (X) = -(L^)- 1 = -(L + i* Y (P° , Vl -P*,M~\ 
and we have the Krein-type resolvent formulas: 

(25) {A- A)- 1 - (A, - A)- 1 = K^ x {L^)-\K'~\y = -K^ x M L {\){K'~\y . 

In the case of a Neumann-type boundary condition as in (22), L x = C — P^ (uVl 
on D{L X ) = D(L). 

8. Applications of pseudodifferential methods I: 
Conditions for lower boundedness 

The formulas we have shown so far use the terminology of V'dbo's mainly as a 
way to indicate what the ingredients in certain operator compositions are. The next 
question to consider is how properties of A are reflected in properties of L. Part of 
the analysis can be carried out with methods of functional analysis, but there also 
exist problems that are solved most efficiently by involving deeper pseudodifferential 
principles. 

An example of how functional analytic and pseudodifferential methods are useful 
together, is the question of lower boundedness inequalities. 

We here restrict the attention to the symmetric set-up where A is formally self- 
adjoint (so A max — Aj^ in ) and Ay is selfadjoint; methods for extending the results 
to nonsymmetric set-ups are found in [G74]. We assume that A 1 has positive lower 
bound m(A 7 ) (cf. (11)). 

In the symmetric case, the general extensions A can of course be nonsymmetric 
(since A max is so). Let us speak of the "selfadjoint case" when only selfadjoint A 's 
are considered. 

In the following, we assume throughout that A corresponds to T: V — > W as 
in Theorem 5.1, and to L: X — > Y* as in Theorem 6.2. The lower boundedness 
problem is the problem of how information on lower bounds on A is related to 
similar information on T or L. The following was shown in [G70] (together with 
studies of coerciveness estimates): 

Theorem 8.1. In the symmetric set-up with A 1 selfadjoint positive, let A corre- 
spond to T : V —¥ W as in Theorem 5.1. 

1° IfVcW and T has a lower bound m(T) satisfying m(T) > —m(A~ ( ), then 
m{A) > m(T)m(A 7 )/(m(T) + m(A y )). 

2° Assume that A 1 is the Friedrichs extension o/^4 m ; n . If m(A) > — oo, then 
V C W andm(T) > m(A). 

In the selfadjoint case these rules go back to Birman [B56], preceded by sesquilin- 
ear form results by Krein [K47] . 

The properties of T are easily translated to similar properties of L using the 
homeomorphism (6); here when X C Y, we set 

(26) m_i(L) = inf{Re(L<p,<p) Y ',Y I V & D(L), \\<p\\_i = 1}, 
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for some choice of norm ||<p||_i on H 2 (£) . (One can let 72 be an isomctry, to 
carry numerical information over between T and L. Sometimes qualitative objects 
such as the sign of m_i(L) are sufficiently interesting.) To take A 1 as the Friedrichs 
extension of A m { n means that it is taken as the Dirichlet realization. 

In 1° we see that the statement "m(T) > —00 ==> m(A) > —00" holds under 
the additional assumption that m(T) > — m{A~ l ); there is a nontrivial question of 
when that assumption can be removed. In [G74] it was shown that when A 1 is the 
Friedrichs extension and A~ x is a corn-pact operator, then m(T) > —00 docs imply 

m(A) > —00. This same result was also announced by Mikhailets and Gorbachuk 
in [GM76] for the selfadjoint case. 

In the application to boundary value problems, we therefore have from this early 
result that lower boundedness of A and L hold simultaneously when we consider 
problems on bounded domains f2, for then A^ 1 is indeed compact. 

For unbounded domains, the question has, to our knowledge, remained unsolved 
up until recently. The question is closely connected with the comparison of T with 
T A as in Theorem 5.3. Indeed, as shown in [G74]: 

Proposition 8.2. Let Gy W be as defined in Theorem 5.3. The property 

(27) m{G zz ) -» 00 for fi ->• -00, 
is necessary and sufficient for the validity of 

(28) m(T) > -00 => m(A) > -00 
for general closed A G M. 

The question was also studied later by Derkach and Malamud [DM91] who 
worked out an analysis that generalizes Proposition 8.2 and gives further condi- 
tions for the validity of the conclusion from m{T) to m{A). However this did not 
capture elliptic problems on unbounded domains (n > 2). 

Because of the recent interest in the analysis of extensions, we have considered the 
problem again, and found a solution in [Gila] for exterior domains (complements 
in M" of bounded domains). 

Theorem 8.3. Let ft be the complement of a smooth bounded set 0_ in R n , and 

let A be symmetric and uniformly strongly elliptic on fl with coefficients in C£°(Q), 
and with a positive lower bound for A 7 . In the application of the extension theory 
to this situation, (27) holds, and hence also (28). 

Here C^°(Cl) stands for the C°°-functions that are bounded with bounded deriva- 
tives. 

The proof relies on the "translation" of abstract operators A to concrete oper- 
ators defined by boundary conditions. Indeed, it turns out that the lower bound 
of G zz behaves like the lower bound m_i(Q^) (cf. (26)) of = P° o Vi - P^. 
Then the deep part of the proof lies in setting the operator in relation to the 
analogous operator for the interior domain f2_, Qt, which does have the desired 
property, in view of our knowledge of problems on bounded domains. The point is 
to show that — Qt)ip, <p)\ is bounded by cH^Hf, / E ) uniformly for fi — > —00, 

so that addition of — Q^L to Qt does not violate the growth of the lower bound. 
This goes by a delicate application of the Vdbo calculus. Details are in [Gila]. 

Let us mention that there is a considerably easier result that holds regardless of 
boundedness of dft and only requires some uniformity in the estimates of coefficients 
of the operators, namely preservation of coerciveness inequalities (Garding inequal- 
ities). We here assume that O is a subset of R™ with smooth boundary, admissible 
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as defined in [G96] (besides bounded domains, this allows exterior domains, per- 
turbed halfspaces and other cases that can be covered with a finite system of local 
coordinates of a relatively simple kind). Moreover, we assume that A is uniformly 
strongly elliptic on ft with coefficients in C£°(f2). The result is that then A satisfies 
a Garding inequality (with c > 0, k G R) 

Re(Au,w) > c\\u\\l - k\\u\\l, u e D(A), 

if and only if X C Y and L satisfies an inequality 

(29) Re(L^,^)> c '||^||^ (s) -/c'||^|| 2 ff _ i(s) , <p e D(L). 

In the case of differential or pseudodifferential Neumann-type boundary condi- 
tions, the inequality (29) for L holds precisely when the pseudodiffcrcntial operator 
it acts like, is strongly elliptic. Details for bounded sets are in [G70] for realizations 
of 2m-order operators; the extension to unbounded sets is shown in [Gila] — the 
argumentation just involves standard trace theorems and interpolation inequalities. 

9. Applications of pseudodifferential methods II: 
Spectral asymptotics 

In the symmetric set-up, when is bounded, the eigenvalues of the selfadjoint 
operator A 7 form a sequence Xj going to oo on R. Already in 1912, Weyl [W12] 
showed the famous estimate for A = — A, for n = 2, 3, m — 2: 

(30) A,(A 7 ) - c 0] m/n is o{j m ' n ) for j ^ oo, 

where Co is a constant defined from the volume of f2; the eigenvalues are repeated 
according to multiplicities. Equivalently, the counting function N(t; A~ ( ) (counting 
the number of eigenvalues in [0,i]), and the eigenvalues fij of the inverse A" 1 , 
satisfy 

N(t; A 7 ) - c A t n l m is o{t n ' m ) for t -> oo, 
^(A- 1 ) c™ /n r m/n is o{r m/n ) for j oo, 

where 

ca = dxd£ 

Jxen, a°(x,£)<l 

(and Co = c A m ^ n ). The estimates have been shown for general n and sharpened 
since then, with more precision on the remainder, and the validity has been ex- 
tended to general elliptic operators A and boundary conditions, and to elliptic 
pseudodifferential operators P on compact manifolds. These improvements have a 
long history that we shall not try to account for here in detail; they are interest- 
ing not only because of the results but even more because of the refined theories 
that were invented in connection with the proofs (for example: Fourier integral 
operators). See e.g. Hormandcr [H68], [H71]. Briining [B74], Seeley [Se78], Ivrii 
[182, 184, 191], Safarov and Vassiliev [SV97]. 

We shall here be concerned with a slightly different question, namely the spectral 
behavior of the difference between the resolvents of two realizations of A. 

Birman showed in [B62] for second-order symmetric uniformly strongly elliptic 
operators A that the singular numbers Sj(B) = [ij(B*B)z of the compact operator 
B = A^ 1 — Az 1 satisfy an upper estimate: 

(31) SjiA- 1 - AZ 1 ) < Cr 2/(n_1) , for all j, 

when A is a selfadjoint realization of A defined by a Neumann or Robin condition. 
In other words, A -1 — Az 1 belongs &{{ n -\)/2), where &( p ) is the space of compact 
operators B for which sj(B) is 0(J~») (often called a weak Schatten class). It is 
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particularly interesting that Birman showed this not just for interior, but also for 
exterior domains, and under low smoothness assumptions. 
We note in passing that the estimate (31) implies that 

(32) A' 1 - A' 1 e C p for p > (n - l)/2, 

where C p is the p-th Schatten class (consisting of the operators B such that 
(sj(B))j e fi 6 £ p (N)); but this is less informative that (31). 

One of the fundamental ingredients in these studies is embedding properties, 
more precisely the knowledge that an operator B that is continuous from L 2 (£l) 
to H s (fl) for some s > (fi bounded smooth C W 1 ) is compact in £ 2 (£!) an d 
belongs to &r n / s \. However, this alone only gives upper bounds on the behavior of 
singular numbers. To get Weyl-type limit properties one must know more about 
the differential or pseudodifferential structure of the operators. 

The estimate (31) was sharpened to a Weyl-type asymptotic estimate in joint 
works of Birman and Solomiak 1978-80, see in particular [BS80] which showed a 
general principle for the spectrum of a ratio of two quadratic forms, implying that 

(33) s^A- 1 - A-'lf^-V -> c, for j -> oo, 

for interior and exterior smooth domains. 

Prior to this, a far-reaching result had been shown in [G74], Section 8: We con- 
sider a symmetric, strongly elliptic 2m-order operator acting in an iV-dimensional 
vector bundle E over a smooth compact Riemanninan manifold VL with boundary 
E, assuming that the Dirichlet realization A 1 is invertible. Let Ab be a selfadjoint 
invertible realization defined by a normal boundary condition 

(34) J2 B oklkU = 0, j = 0, 1, . . . , 2m - 1, 

k<j 

where the Bjk are differential operators of order j — k from to given vector 
bundles Fj over S (with dim .Fj > 0); normality of the boundary condition means 
that the Bjj are surjective morphisms. (That A acts in the vectorbundle E means 
that it is locally (N x 7V)-matrix-formed. In the scalar case, N = 1 and the Fj are 
0- or 1-dimensional, with Bjj an invertible function when dim Fj = 1. Ellipticity 
of the boundary condition requires in particular that ^\ dim .Fj = mN.) Denote 
®j> m Fj — F 1 . 

Theorem 9.1. Let T: V — > V be the operator corresponding to Ab by Theorem 
5.1. There exists an isometry J : L 2 (E,F 1 ) ~> V with inverse J^ 1 = J* (in the 
Tpdbo calculus), such that 

Ti = J*TJ 

acts like an elliptic invertible ipdo 7~ in F 1 of order 2m, and D(7i) — 
{<fi E L 2 (T,, F 1 ) | T(fi G L 2 (Z,F 1 )} = H 2m (Y,, F 1 ). Here Ti has the same spec- 
trum as T, and its eigenvectors are mapped to the corresponding eigenvectors ofT\ 
by the isometry J . Moreover, 

(35) A B X - A^ 1 = \ v JTi X J* pr y , 
whereby the positive resp. negative eigenvalues satisfy 

lifiAz 1 -A- 1 )= t ,f(T 1 - 1 ) for all j. 
It follows that with constants determined from the principal symbols, 
N !± (t; A B X - A- 1 ) = c^"- 1 )^ 2 ™) + o( t (n-i)/(2m)) for t 

flfiA-i-A- 1 ) = (C±) 2 W(»-l) r 2m/(n-l) +0(r (2m + l)/(n-l) ); 

here N r± (t; S) indicates the number of positive, resp. negative eigenvalues of S 
outside the interval } — 1/t, l/t[ . 
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The two statements in (36) are equivalent, cf. e.g. [G96], Lemma A. 5. They follow 
from Hormander [H68] for elliptic ^do's, when the principal symbol eigenvalues of 
T are simple; this restriction is removed by results of Ivrii [182] . See also Theorem 
10.1 below. 

We have recently checked that the proof extends to exterior domains, for uni- 
formly strongly elliptic systems. 

The fine estimates with remainders depend on the ellipticity of the ^>do. For 
simple Weyl-type estimates, Birman and Solomiak in [BS77] removed the ellipticity 
hypothesis, showing that 

Sj {P)j k/n ->c(p°) for j -> oo 

holds for any classical ^do P of order — k < on a compact manifold of dimension n. 
They even allowed a certain nonsmoothness of the homogeneous principal symbol, 
both in x and £, needing only a little more than continuity. In the elliptic case, there 
are recent works of Ivrii dealing with remainder estimates under weak smoothness 
hypotheses. 

In [G84] we made an effort to increase the accessibility of the ^dbo calculus 
by publishing an introduction to it with several improvements, and showing as a 
main result that any singular Green operator of negative order and class has a 
Weyl-type spectral asymptotics formula: 

Theorem 9.2. Let G be a classical singular Green operator of order —k < and 
class on an n- dimensional compact manifold with boundary. It has a spectral 
asymptotics behavior 

(37) Sj (G)j fe /("- 1 ) -+ c(g°) for j oo, 

where c(g ) is a constant defined from on the principal symbol of G. 

This was moreover used to show asymptotic formulas generalizing (33), both 
for interior and exterior domains, followed up in another study [G84a] including 
also the dependence on a spectral parameter A. Indeed, we have as an immediate 
corollary of Theorem 9.2, also for nonselfadjoint cases: 

Corollary 9.3. Let A be elliptic of order 2m with invertible Dirichlet realization 
and let As be an invertible realization defined by a normal elliptic boundary condi- 
tion. For any positive integer N, A^ N — A~ N is a singular Green operator of order 
-2mJV and class 0, and hence satisfies 

(38) s 3 {A B N A - N ).f mN ^-V CN for j -+ oo, 

for a constant cn defined from the principal symbols. 

Also exterior domains are considered in [G84], [G84a], where (38) is shown for 
realizations of second-order operators and their iterates. The results apply of course 
to resolvents by replacement of A by A — A; the A-dependence is studied in [G84a]. 
It is seen that the tpdbo theory provides a forceful tool for such questions, and we 
strongly recommend its use. 

10. New spectral results 

Spectral estimates of resolvent differences have been taken up in recent papers by 
Behrndt et al. [AB09, BLLLP10, BLL11] for second-order operators and Malamud 
[M10] for 2m-order operators with normal boundary conditions, based on boundary 
triples methods. Here Schatten class and weak Schatten class estimates are shown, 
relying on such estimates for Sobolev space embeddings. 

We have returned to the subject in [Gil] where we, besides showing new results 
on perturbations of essential spectra, have reformulated and extended results in 
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[G84] on estimates like (38), including general differences and exterior domains. 
The central ingredient is the estimate (37) for singular Green operators, plus the 
fact that s.g.o.s give their essential contribution in a small neighborhood of the 
boundary, also for exterior domains, allowing cutoffs eliminating infinity. 

An inspection of the results of [G74] shows that the spectral estimates in Theorem 
9.1 can be further sharpened by use of results of Ivrii [182]: 

Theorem 10.1. In the setting of Theorem 9.1, assume in addition that the prin- 
cipal symbol of T satifi.es Ivrii 's conditions (H± ) from [182] ( the bicharateristics 
through points o/T*(X) \0 are nonperiodic except for a set of measure zero). Then 
there are constants Cf such that 

(39) N^it; A B X - A- 1 ) = c± t ("-i)/(2»0 + ^(M)/(H + (£(™- 2 )/( 2 ™)). 
The proof is a direct application of [182] Th. 0.2 to T ■ 

The formula (35) is a special type of Krein resolvent formula with isometrics, 
valid for selfadjoint realizations, but the analysis in [G74] also implies Krein for- 
mulas in the nonselfadjoint cases. Namely, Th. 6.4 there shows how T: V —¥ W 
is represented by a realization C\ of a ip&o C acting between vector bundles over 
S, and here Ab is elliptic if and only if C is elliptic (Cor. 6.10). In the invertible 
elliptic case, formula (20) then takes the form 

A B l - A- 1 = Wlv^^i^wT 1 Ww, 

with £i acting like -(B w + P n P 7:X )$ (notation explained in [G74]); the right- 
hand side is a composition of a Poisson operator, an elliptic V'do and the adjoint of 
a Poisson operator, all of mixed order. Its s-numbers can be studied by reduction 
to an elliptic ?/>do over E, where Ivrii's sharp results can be used. 

Let us just demonstrate this for second-order operators, for the formula (23), 
with L- 1 elliptic of order -1: Denote A- 1 - A^ 1 = S. Then 

Sj (Sf = Sj (K 7 L-\K' 7 rf = ^K^iK^K^L-'YK;) 
= H j (L-\K'J*K' 7 (L- 1 rK;KJ, 

where we used the general rule Hj(BiPt2) = /Uj(i?2-Bi). Both operators P\ — K*K 7 

and P[ = (K'^YKly are selfadjoint positive elliptic V'do's of order —1 (cf. e.g. [Glib], 

i 

proof of Th. 4.4). Let P 2 = P^ , then we continue the calculations as follows: 

Sj {Sf = ^(L-ipKL-'yP^ = ^■(P 2 L- 1 P 1 '(L- 1 )*P 2 ) = ^(P 3 ), 

where P3 = P2i _1 P 1 '(i _1 )*p2 is a selfadjoint positive elliptic -0do on S of order 
—4. Applying Ivrii's theorem to P 3 _1 , we conclude: 

Theorem 10.2. For the operator considered in Theorem 6.3, the s-numbers satisfy 

(40) N'(t; A' 1 - A' 1 ) = Ct^-W 2 + 0(i ( ™- 2)/2 ). 

Moreover, if the principal symbol of P 3 _1 satifies Ivrii's condition from [182] (the 
bicharateristics through points o/T*(£)\0 are nonperiodic except for a set of mea- 
sure zero), there is a constant C\ such that 

(41) N'(t; A- 1 - A- 1 ) = Ct^l 2 + C^"" 2 )/ 2 + o(^"- 2 )/ 2 ). 

Sharpened asymptotic formulas can also be obtained for differences between re- 
solvents of two realizations that both differ from the Dirichlet realization, by use of 
the analysis in [G68] with a general invertible realization Ap as reference operator. 

To give another example of applications of the V'dbo theory, the following result 
is found straightforwardly as a consequence of [G84]: 
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Theorem 10.3. Let Ab and Ab> be elliptic invertible realizations of A such that 
B and B' map into the same bundles and have the same principal part. Then 
Ag 1 — Ag} is a singular Green operator of order —2m — 1 (since its principal part 
is zero), and hence, by (37), 

s^Ag 1 - J 4", 1 )j( 2 '"+ 1 )/("- 1 ) -> c for j -> oo. 

The singular Green operator will be of a still lower order —2m — r if the first 
r > 1 terms in the symbols of B and B' coincide. 

Example 10.4. As a special case, we can compare two Robin conditions for a 
second-order operator A: 

A\ defined by v\u = bijou, 

A2 defined by v\u = &270U; 

&i,&2 £ C°°(E). When regarded from the point of view of Theorem 10.3, these 
are normal boundary conditions v\u — B^qu — 0, where B\ = bi and B 2 = b 2 
considered as first-order operators have principal part 0, so the boundary operators 
have the same principal part. Then the s.g.o. A^ 1 — A^ 1 is of order —3, and by 
(37), 

(42) Sj (A^ - A^f'^-V -> c for j -> 00. 

In [BLLLP10], Berndt, Langer, Lobanov, Lotoreichik and Popov prove upper 
estimates for this difference in the case A = — A— A, namely Schatten class estimates 
of s^A^ 1 - A^ 1 ) as in (32) with (n - l)/2 replaced by (n - l)/3. 

In case the bi are C°°, the result is covered by (42) as explained above. However, 
the bi in [BLLLP10] are allowed to be nonsmooth, namely to be in L^iT,), which 
goes outside the range covered by the smooth -0dbo theory. 

This led us to investigate how far we could push the proof of asymptotic estimates 
(42) to make them valid for nonsmooth choices of bi. The outcome is published 
in [Glib], where it is shown that (42) holds for symmetric second-order strongly 
elliptic operators on smooth domains, when b\ and b 2 are piecewise C £ on E, having 
jumps at a smooth hypersurface. 

The Schatten class estimates have been followed up by Berndt, Langer and Lo- 
toreichik in a study of sclfadjoint realizations [BLL11]. 

Unsolved questions of asymptotic estimates lie primarily in the range of sit- 
uations with limited smoothness. Resolvent formulas have been studied in such 
general cases, [PR09] and [G08] for C M -domains, [GM08, GM11] for Lipschitz and 
quasi-convex domains, [AGW11] for a class of domains containing C 3//2+£ , with a 
nonsmooth generalization of ^dbo's. To our knowledge, spectral asymptotic esti- 
mates have not yet been worked out for such resolvent differences. Some upper 
estimates are in sclfadjoint cases known from Birman [B62]. 

A problem with a different flavor is the case of a mixed boundary condition, such 
as prescribing for —A the Dirichlet condition on a part E_ of the boundary and 
a Neumann- type condition on the other part E + . Here there is a jump in the 
order of the boundary condition. The domain of the realization is contained in 
H 2 -e (f2) only for e > 0, so it is not covered by boundary triples methods requiring 
the domain to be in i?5(Q). Spectral upper estimates are known from [B62]. A 
spectral asymptotic estimate was obtained recently in [Gllc], based on somewhat 
intricate applications of results on nonstandard pscudodiffcrcntial operators. 

There are also other questions that can benefit from pseudodifferential meth- 
ods, for example the study of spectral asymptotics of the nonelliptic Krem-like 
extensions, cf. [Gila]. 
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